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ABSTRACT 

The spatial origin and detectability of rotational H 2 emission lines from Herbig Ae type protoplanetary disks beyond 70yum is 
discussed. We use the recently developed disk code ProDiMo to calculate the thermo-chemical structure of a Herbig Ae type disk 
and apply the non-LTE line radiative transfer code Ratran to predict water line profiles and intensity maps. The model shows three 
spatially distinct regions in the disk where water concentrations are high, related to different chemical pathways to form the water: 
(1) a big water reservoir in the deep midplane behind the inner rim, (2) a belt of cold water around the distant icy midplane beyond 
the "snowline" r i, 20 AU, and (3) a layer of irradiated hot water at high altitudes z/r = 0.1 ...0.3, extending from about 1 AU to 
30 AU, where the kinetic gas temperature ranges from 200 K to 1500 K. Although region 3 contains only little amounts of water 
vapour (~ 10~ 4 Msai-thX it is this warm layer that is almost entirely responsible for the rotational water emission lines as observable 
with Herschel. Only one ortho and two para H 2 lines with the lowest excitation energies < 100 K are found to originate partly 
from region 2. We conclude that observations of rotational water lines from Herbig Ae disks probe first and foremost the conditions 
in region 3, where water is predominantly formed via neutral-neutral reactions and the gas is thermally decoupled from the dust 
Tgas > Tdust- The observation of rotational water lines does not allow for a determination of the snowline, because the snowline 
truncates the radial extension of region 1, whereas the lines originate from the region 3. Different line transfer approximations (LTE, 
escape probability, Monte Carlo) are discussed. A non-LTE treatment is required in most cases, and the results obtained with the 
escape probability method are found to underestimate the Monte Carlo results by 2% -45%. 
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1. Introduction 

Water is one of the most important species in planet formation, 
disk evolution and for the origin of life. In protoplanetary disks, 
water can be abundant either in the gas phase or as solid ice, 
owing to a high sublimation temperature. 

Water vapour is predicted to be abundant inside of the snow- 
line within a few AU where densities are high and temper- 
atures are too warm f or ice formation (Agu ndez et all [2008; 
Glassgold et al.L I2009I) . Ind eed, observa t ions i n the near- and 
mid-infrared with Spitzer (ISalvk et all 120081 ICarr & Naiital 
200& lNaiita et all 120091; lEisner et al.ll2009l) reveal the presence 
of water vapour and other simple organic molecules (OH, C2H2, 
HCN, CO2). Simple modelling of the line emissions constrains 
the warm gas (a few 100 K) to be located inside of r ^ 3 -5 AU. 

The core-accretion model for planet formation is only ef- 
ficient in the cold midplane of disks where dust grains are 
covered by water ice. Icy grains are stickier than bare silicat e 
grains and coagulate faster into planetesimals (llda & Linl l2008). 
Liquid water is one of the prerequisites for the emergence of life 
on terrestrial planets. But its origin, either via release of water 
molecules trapped in hydrated rocks dur ing volcanis m or via the 
impact of comets, is still being debated (iNu thi 120081) . 

The lowest rotational lines of water lie in the far IR and 
are only observable by satellites, e. g. the Herschel Space 



Observatory. Other recent works on rotational water lines from 
protoplanetary disks used X-ray models (without UV pho- 
toprocesses) to calculate the thermo-chemical disk structure 
on to p of a pre-described density structure (Meii erink et all 
120081) and applied a mu l ti-zone escape probability method 
(Poelman & Spaansi I2005L |2006) to compute the line fluxes. 
In this Letter, we use the disk code ProDiMo to compute 
the disk structure, temperature, and water abundance self- 
consistently, and discuss the prospects for detecting rota- 
tional water lines by means of th e Monte Carlo code Ratran 
( Hogerheij de & van der TjkLl2000). 

2. The model 

We used the recently developed disk code ProDiMo to cal- 
culate the thermo-chemical structure of a protoplanetary disk 
around a Herbig Ae type star with parameters listed in Table [T] 
ProDiMo combines frequency-dependent 2D dust-continuum 
radiative transfer, kinetic gas-phase and UV photo-chemistry, 
ice formation, and detailed non-LTE heating & cooling with the 
consistent calculation of the hydrostatic disk structure. ProDiMo 
does not include X-rays at the moment. X-ray to FUV lumi- 
nosity ratios of He rbig Ae stars are often low Lx/^fuv ^ 1 
dKamp et all 12008b . so that we assume that the FUV irradiation 
provides the main energy input for the disk. The model is char- 
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Fig. 1. Concentration of water molecules en 2 o = «h 2 o/ w <h> as 
function of radius r and relative height above the midplane z/r 
for a Herbig Ae disk model. Three different regions with high 
water concentration can be distinguished as described in the 
text. The two white contours show r gas = 200 K and 1500 K. 
The red contour line marks total hydrogen nuclei particle den- 
sity n <H > = 10 10 - 7 cnr 3 . 



Table 1. Herbig Ae type disk model parameter. 



quantity 


symbol 


value 


stellar mass 




2.2 M 


effective temperature 




8600 K 


stellar luminosity 


L+ 


32 L 


disk mass 


M disk 


0.01 M Q 


inner disk radius 


«m 


0.5AU™ 


outer disk radius 


^out 


150 AU 


radial column density power index 


e 


1.0 


dust-to-gas mass ratio 


Pdlp 


0.01 


minimum dust particle radius 


^min 


0.1 fim 


maximum dust particle radius 


^max 


200 fim 


dust size distribution power index 
dust material mass density <2) 


^pow 


3.5 


Pgr 


2.5 g cm -3 


strength of incident ISM UV 
abundance of PAHs relative to ISM 


^ISM 

/pah 


1 

0.1 



(1) : soft inner edge applied, see Woitk e et al. (200 9) 

(2) : dust optical constants from lDraine" & Led fl 984) 



acterised by a high degree of consistency between the various 
physical, chemical, and radiative processes, where the mutual 
feedbacks ar e solve d by global iterations. For more details see 
IWoitke et alj d2009l) . henceforth called Paper I. Recent updates 
include an improved treatment of UV photorates by detailed 
cross sections in the calculated radiation field (see iKamp et aU 
120091) . 

The model results in a flared disk structure with a puffed- 
up inner rim and a vertically extended hot atomic layer above 
z/r ^0.15-0.4 from the inner rim to aboutr = 20AU, similar to 
Fig. 9 in Paper I (l.h.s.), where r gas > 1000 K due to the stellar 
UV-irradiation. The major difference between the T Tauri type 
disk discussed in Paper I and the Herbig Ae disk discussed here 
is that in Herbig disks, the star is much more luminous, so the 
dust is warmer in the midplane (here 7d ust > 100K inward of 
r — 20 AU), which prevents water ice formation. 

Water molecules generally form in deeper layers, and the re- 
sulting water concentration in these layers is depicted in Fig.Q] 



Table 2. Characteristics of water regions in Herbig Ae disk 
model. 



region 


r[AU] 


z/r 


Mn,o [M Earth ] 


(T gas ) [K] 


(W [K] 


1 


0.7 - 20 


<0.1 


0.31 


240 


240 


2 


20-150 


<0.1 


3.2 x 10~ 5 


38 


37 


3 


1-30 


0.1-0.3 


1.1 x 10~ 4 


410 


150 



The vertical H2O column densities in this model are found to 
be 10 22 cnr 2 at 1 AU, still 10 19 cnr 2 at 10 AU, but then quickly 
dropping below 10 15 cirT 2 at 30 AU and beyond. 

3. Chemical pathways to water 

The formation of H2O follows different chemical pathways in 
the three different regions shown in Fig.Q] Two of these regions 
(1 and 3) have been previously ident ified in vertical slab m odels 
for X-ray irradiated T Tauri disks by Glassgold et al. (2009). 

1) The big inner water reservoir. The deep midplane regions 
from just behind the inner rim to a distance of about 10AU in 
this model host the majority of the water in the disk, see Table|2] 
This region is almost completely shielded from the stellar and in- 
terstellar radiation {Ay > 10), is still too warm for water to freeze 
out (7dust > 120 K), and is characterised by large particle densi- 
ties > 10 11 cirT 3 and thermal equilibrium between gas and dust 
(Tgas ~ Tdust)- The important chemical reactions in this region are 
in detailed balance with their direct reverse reaction; i. e. the gas- 
phase chemistry is close to thermochemical equilibrium. Since 
H2O is the thermodynamically most stable oxygen carrier un- 
der these conditions (after CO), the water concentration is given 
approximately by the element abundance of oxygen minus the 
fraction thereof bound in CO. 

2) The distant water belt. Region 2 at r « 20 - 150 AU and 
z/r ^0.05 is characterised by particle densities ~ 10 9 -10 10 cm -3 , 
temperatures r gas « 7d ust « 30- 100K, and UV-strengths y ( see 
Eq. (41) in Paper I) below 500. These conditions allow water to 
freeze out, and there is a tight equilibrium between water ad- 
sorption and photodesorption H2O «-» H20 # , where H20* desig- 
nates water ice. The formation of gaseous water in this region is 
mainly controlled by the following two photoreactions 

H 2 # + hv — > H 2 (1) 
H 2 + hv — > OH + H . (2) 

Therefore, the strength of the UV field x controls the water 
vapour concentratiorQ. At r — 80 AU, x drops from ~ 10 3 to 
virtually between z/r ^0.15 and 0.02. The UV photons mostly 
come from the star and are scattered around the optically thick 
core of the disk. Water is photo-dissociated for^ that is too high, 
and freezes out for too low x- The result is a thin belt of cold 
water vapour with c oncentrations in excess of 10~ 9 around the 
distant icy midplane. Glassg oTd et alj d2009h have not found re- 
gion 2 in their T Tauri models, probably because their model did 
not include the UV photoprocesses. 

3) The hot water layer. At distances r w 1 - 30 AU and relative 
heights z/r 5= 0. 1-0.3, the model is featured by an additional layer 
of warm water-rich gas that is thermally decoupled from the dust 
2gas > ?dust- The particle density is about ~ 10 s - 10 10 crrT 3 in 



1 Our analysis is restricted to the case of kinetic chemical equilib- 
rium. Thus, we cannot discuss the history of water ice formation with 
our model. Water ice may furthermore b e photodesorbed directly into 
OH, see e.g. lAndersson & van Di shoeck (20081). 




Fig. 2. Monte Carlo simulations of three ortho H2O lines with increasing excitation energy for a distance of 140 pc and inclination 
30° (infinite resolution). The inserted maps show the integrated line intensity I(x,y) = J I v (x,y) - Iy° nt (x,y) dv as function of sky 
position {x,y) with 7™ nt (x,y) = I v {x,y) far away from line centre. An angular distance of 0.3" in the maps corresponds to a spatial 
distance of 42 AU. The first row shows the results for the full thermo-chemical input model. The second row shows the results with 
the water molecules being removed from region 3 (see Fig.[TJ. The low-excitation line on the l.h.s. is depicted on a 3x larger angular 
scale. The scaling of the spectral flux F v is partly enlarged by a factor x5 in the lower panels as indicated. The colour-scale for the 
intensity maps is the same in each column. 



this layer, the dust temperature 100-200K, the gas temperature 
200- 1500K, and the UV radiation field strength X ~ 10 4 - 10 7 . 

Above region 3, the high gas temperatures in combination 
with the direct UV irradiation from the star efficiently destroys 
all OH and H2O. In region 3, the medium is shielded from the di- 
rect stellar irradiation by the puffed-up inner rim. In the shadow 
of the inner rim, x drops quickly by about two orders of mag- 
nitude (the remaining UV photons are scattered stellar photons), 
and water forms via the following chain of surface and neutral- 
neutral reactions 

H + H + dust H 2 + dust (3) 
H 2 + O — > OH + H (4) 
OH + H 2 H z O + H , (5) 

counterbalanced by the photo-dissociation reactions OH + hv — » 
O + H and H 2 + hv — > OH + H. The reason for the mini- 
mum of the H 2 concentration between regions 3 and 1 lies in 
the high activation energies required for reactions © and ©. 
These neutral-neutral reactions become inefficient when the gas 
temperature drops below r gas ^ 200 K. The radial extension of 
the hot water layer is r estricted for the same r eason. Since UV 
chemistry is lacking in dGlassgold et all 120091) . their water for- 
mation is counterbalanced by charge transfer reactions with H+, 
which also destroys OH and H 2 0. 



4. Spectral appearance of rotational water lines 

Having calculated the density structure, the molecular abun- 
dances, the dust and gas temperatures and the continuous radia- 
tion field in the disk, we performed axisymmetric non-LTE line 
transfer calculations for selected rotational water lines as sum- 
marised in Table [3] The non-LTE input data for ortho (para) H 2 Q 
is taken from the Leiden Lambda database (Schoi er et all 120051) . 
which includes 45 (45) levels, 158 (157) lines, and 990 (990) col- 
lisional transitions with H 2 . A scaled version of the last data is 
also applied to collisions with H, which is essential as region 3 
is partly H 2 -poor and H-rich. The velocity field is assumed to 
be Keplerian. We add a turbulent line width of 0.15km/s to the 
thermal line width throughout the disk. The ratio between para 
and ortho H 2 is assumed to be as in LTE. 

We used three different methods of increasing complex- 
ity to calculate the water population numbers: local thermal 
equilibrium (LTE), a simple escape probability method (ES, 
see Sect. 6.1 of Paper I) a nd a modified version of the 2D 
Mon te Carlo cod e Ratran dHogerheiide & van derTakL |2000), 
see dKamp et all [2009) for modifications. The LTE and ES 
methods used the full 150x150 ProDiMo output directly as 
thermo-chemical input model. The grid size of the MC model 
needed to be somewhat reduced for practical reasons. We de- 
cided to run 80x80 MC models, which need about 7 x 10 5 pho- 
ton packages to converge to a signal/noise ratio better than 5 for 
the worst population number in the worst cell, which takes about 
13 CPU hours on a 2.66 GHz Linux machine. 
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To investigate the role of the hot water layer (region 3) for 
the spectral appearance of the rotational lines, we calculated two 
sets for each model. The first set includes the full chemical input 
model. In the second set, we artificially put the water abundance 
in region 3 to zero (if n^j) < 10 10 7 cirT 3 and T gas > 200 K). The 
results of the two sets of MC models are compared in Fig. [2] 

We generally observe double-peaked line profiles typical for 
rotating gas in emission. In the case of the full input model, the 
peak separation generally measures the radial extension of re- 
gion 3 (Kepler velocity is 8.8km/s at 25 AU in this model, in- 
clined to 30° gives 4.4km/s). However, the peak separation of 
the lowest 0-H2O excitation line (538.29 fim) and the two lowest 
P-H2O excitation lines (269.27 jum and 303.46/im, see Table [3} 
correspond to the full radial extent of the model, 150 AU. 

The truncated model generally results in 10-20x smaller 
line fluxes with a broader, often unclear profile. If the water in 
region 3 is missing, the lines originate mainly from region 1, 
which is optically thick even in the continuum. Since the gas 
is in thermal balance with the dust in region 1 (T gas « T<i m ), the 
lines do not go much into emission in region 1 . Again, the lowest 
three excitation lines behave differently, and region 2 contributes 
by ~ 30% for these lines. There is furthermore one intermittent 
case (0-H2O 179.5 //m with E u = 114K) where the truncated 
model reveals a small contribution of the extended region 2 with 
the character of low-excitation lines (middle column in Fig. [2j. 

All rotational water lines of the Herbig Ae disk discussed in 
this letter are above the lcr detection limit of the PACS spec- 
trometer (as 1-5x1 0" 18 W/m 2 , depending on A). The strongest 
water lines at 78.74 fim and 89.99 fim are above the 5cr detection 
limit. However, the lines sit on a strong dust continuum, which 
possibly complicates the detection by Herschel. 

5. LTE vs. escape probability vs. Monte Carlo 

Figure [3] compares the results obtained by three different 
line transfer methods for the high-excitation para-H20 line at 
89.99 fim. The results for the other lines are listed in Table[3] We 
consider the deviations from the results of the most advanced 
method (MC) as a measure of the quality of the other methods. 
The LTE predictions are generally too high, by up to a factor of 
3.5, although continuum flux, line width, and peak separation are 
similar. Since the densities in region 3 are lower than the criti- 
cal density ~ 10 10 cirT 3 , the levels tend to depopulate radiatively, 
which explains the overpredictions by LTE. Deviations between 
ES and MC are between 2% and 45%, and increase with excita- 
tion energy E u . Our ES method tends to underestimate the line 
fluxes in general. The levels in region 3 are pumped by line ra- 
diation from distant regions in LTE that have larger line source 
functions. This effect is difficult, if not impossible, to be properly 
account for in the ES approximation. 

6. Conclusions 

The rotational water lines from Herbig Ae disks beyond IQfim 
originate predominantly from a warm molecular layer at relative 
altitudes z/r:«0.1-0.3 where H2O is formed via neutral-neutral 
reactions in a thermally decoupled gas (T gas > rd ust ). The more 
distant cold water around the icy midplane, where the ice is pho- 
todesorbed, contributes only to the lowest excitation lines. The 
peak separation of all other lines measures the radial extension of 
the warm molecular layer, which is about 40 AU in the discussed 
model. In contrast, the vast majority of water vapour in the disk 
is situated in the deep midplane, extending from just behind the 



Table 3. Properties of rotational water lines, and calculated line 
fluxes Fi = J(F V - Fl° n[ ) dv for different line transfer methods. 
Listed line fluxes are in units 10~ 18 W/m 2 for distance 140pc 
and inclination 30°. * indicates that this line will be observed by 
the Herschel OT Key Programme Gasps. ES = escape probabil- 
ity, MC = Monte Carlo. 



carrier 


transition 


4L"mJ 


A„,[s-'] 


E„ [K] 


LTE 


ES 


MC 


|ES-MC| 


o-H z O 


lio 


->loi 


538.29 


0.0035 


61 


0.702 


0.531 


0.538 


2% 


o-H 2 


4 2i 


^3,2* 


180.49 


0.0306 


194 


9.27 


2.32 


3.39 


37% 


o-H 2 


2,2 


->lot* 


179.53 


0.0559 


114 


11.9 


7.69 


8.59 


11% 


0-H2O 


3 3 


^2 12 


174.63 


0.0505 


197 


12.7 


7.47 


8.88 


17% 


0-H2O 


2l2 


->lio 


108.07 


0.256 


194 


54.8 


11.4 


17.9 


44% 


0-H2O 


4 23 


^3 I2 * 


78.74 


0.484 


432 


94.5 


24.6 


35.9 


37% 


p-H 2 


2()2 


111 


303.46 


0.0058 


101 


1.27 


1.48 


1.40 


5% 


p-H 2 


111 


->Ooo 


269.27 


0.0184 


53 


2.40 


2.22 


2.09 


6% 


p-H 2 


4,3 


^ 3 22 * 


144.52 


0.0332 


397 


5.53 


4.57 


4.15 


9% 


p-H 2 


3l3 


— * 2 2 


138.53 


0.125 


204 


14.3 


7.76 


8.71 


12% 


p-LLO 




111 


100.98 


0.260 


196 


32.5 


8.66 


10.9 


23% 


p-H 2 


322 


->2n* 


89.99 


0.352 


297 


48.0 


10.2 


14.0 


31% 
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Fig. 3. Comparison between the results of different line transfer 
methods in application to the high-excitation para-H20 line at 
89.99 fim. The full black line redisplays the results of the Monte 
Carlo treatment already shown in Fig. [2] upper row, middle fig- 
ure. 

inner rim outward to the snowline, where water freezes out to 
form water ice. The gas in this massive deep water reservoir is in 
thermal balance (r gas « 7d ust ) with optically thick dust and, there- 
fore, no strong line emissions are produced with respect to the 
continuum from this deep region. Thus, no information about the 
position of the snowline can be deduced from th e rotational wa- 
ter lin es. A similar conclusion was reached by iMeiierink et alj 
d2008l) . The line analysis generally requires a non-LTE treat- 
ment. Our escape probability method is found to underestimate 
the water line fluxes with respect to the more expensive Monte 
Carlo method by about 2% - 45%. 
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